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Abstract 

A general formalism for computing only the rational parts of one- 
loop QCD amplitudes is developed. Starting from the Feynman inte- 
gral representation of the one-loop amplitude, we use tensor reduction 
and recursive relations to compute the rational parts directly. Explicit 
formulas for the rational parts are given for all bubble and triangle 
integrals. Formulas are also given for box integrals up to two-mass- 
hard boxes which are the needed ingredients to compute up to 6-gluon 
QCD amplitudes. We use this method to compute explicitly the ratio- 
nal parts of the 5- and 6-gluon QCD amplitudes in two accompanying 
papers. 
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1 Introduction 



The forthcoming experimental program at CERN's Large Hadron Collider 
(LHC) requires many computations at the next-to-leading order (NLO) or 
one-loop with many particles as final states [1]. However the analytic com- 
putation of the matrix elements is very difficult. Only for special helicity 
configurations [2, 3, 4] or special models [5], some analytic results are known 
for higher point amplitudes. The current state of the art in NLO compu- 
tation is 5-point for QCD processes and 6-point for electroweak processes 
[6, 7, 8]. The recent development in tackling the multi-leg amplitudes by 
semi-numerical/analytic methods shows promise for improving traditional 
capabilities [9, 10, 11, 12, 13, 14, 15, 16]. All helicity configurations for the 
6-gluon amplitude are evaluated for a single-space point [9]. These results 
are used to check the available analytic results [17]. 

Following Witten's twistor string theory [18], the CSW approach [19] and 
the use of maximally-helicity-violating (MHV) 1 vertices [21, 22], there has 
been spectacular progress [23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 
37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47] in the perturbative QCD computations 
in the last two years or so, by using the unitarity cut method of Bern, Dunbar, 
Dixon and Kosower [3, 48, 50] and the spinor-helicity formalism [51, 52] 
(see [53] for a review). In particular, Bedford, Brandhuber, Spence and 
Travaglini [28, 31] applied the MHV vertices to one-loop calculations. Britto, 
Buchbinder, Cachazo, Feng and Mastrolia [38, 39, 40] developed an efficient 
technique for evaluating the rational coefficients in an expansion of the one- 
loop amplitude in terms of scalar box, triangle and bubble integrals (the 
cut-constructible part, see below and Sect. 4). By using their technique, it is 
much easier to calculate the coefficients of box integrals without doing any 
integration. Recently, Britto, Feng and Mastrolia completed the computation 
of the cut-constructible terms for all the 6-gluon helicity amplitudes [40]. 

In order to complete the QCD calculation for the 6-gluon amplitude, 
the remaining challenge is to compute the rational part of the amplitude 
with scalars circulating in the loop, commonly called the N = case in a 
super symmetric decomposition of QCD amplitudes: 

j^qcd A N=4: AA N= ^ chiral -\- A N= ® or scalar (1) 



1 The 2 dimensional origin of the MHV amplitudes in gauge theory was first given in 
[20] 
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The above strategy of splitting the computation of the QCD amplitude into 
various supersymmetric parts plus a scalar part is quite fruitful. By using 
a theorem of Bern, Dunbar, Dixon and Kosower [50], the supersymmetric 
parts are cut-constructible, meaning that these amplitudes can be determined 
completely by using 4-dimensional unitarity. Even for the scalar part which is 
not cut-constructible, we can still split it into two parts: a cut-constructible 
part and a rational part. As we said before, the recent development inspired 
by twistor string theory has lead to very efficient techniques to compute 
the cut-constructible part [38, 39, 40]. To complete the program it is quite 
important to have efficient and powerful methods to compute the rational 
part. 

There are various ways to compute the rational part. The first approach 
[48] is to use the factorization properties by trial and error. This is a quite 
effective method if the final result is simple enough. The difficulty with 
this approach is that we do not know how to automate the method to make 
effective use of the advances in computer industry. The correctness of the ob- 
tained result is almost guaranteed by checking factorizations for all channels. 
For higher point amplitudes, the complexity of the analytic results makes 
this method impractical. 

The second approach uses the unitarity relation. In principle the rational 
part can be constructed by using the D-dimensional unitarity method [49, 
50, 54, 42]. The problem with this approach is that too much information 
is kept and tree amplitudes in D-dimension are even more difficult [55]. In 
fact this approach loses the simplicity of 4-dimensional helicity amplitudes 
as given by the MHV formula. 

The third approach is the bootstrap recursive approach of Bern, Dixon 
and Kosower [34]. This approach is quite promising and powerful. It is a 
streamlined approach of the first one by adding the insights of the recent 
tree-level recursive method of Britto, Cachazo, Feng and Witten [32]. This 
approach has already produced a wealth of general results for special helicity 
configurations, notably the "split-helicity" configurations [34, 46, 17]. It can 
also be used to compute one-loop QCD amplitudes with general helicities as 
outlined in [17]. 

Given the complexity of the results for the cut-constructible part of the 
6-gluon amplitude [40] and its important applications to LHC related ex- 
periments, it is quite worthy and even mandatory to have other methods to 
compute the rational part of the QCD amplitude. In particular one would 
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like to bypass the need of using the cut-constructible part and have an inde- 
pendent method to compute the rational part. Of course, the testing ground 
for any method is a complete computation of the 6-gluon QCD amplitude 
where only partial results for some helicity configurations exist. The present 
status for the marching to one loop 6-gluon QCD amplitude were summarized 
in [55, 56]. For more recent developments, we refer the reader to [57, 58]. 

In this paper we will study the problem of computing the rational parts of 
one-loop amplitudes directly from Feynman integral representations. These 
integrals can be written down directly by drawing all Feynman diagrams and 
by using the Feynman rules. With the present technology these can be done 
quite effectively by using the various packages like GRACE [62], FeynArts 
[63] and Qgraf [64] et. al. (See [65, 62] for reviews.) Fortunately these 
powerful methods are not needed to compute up to the 6-gluon amplitude. 
For computing higher point amplitudes or 6-parton amplitudes they may be 
a necessity. 

It is easy to imagine that the rational part is already contained in the 
integral representation of the amplitude. If one could obtain the complete 
rational coefficients by doing tensor reduction to scalar box, triangle and 
bubble integrals, one can simply get the rational part by making an expansion 
with the dimension D around 4 (D = 4 — 2e is the parameter of dimensional 
regularization). This is extraordinarily difficult because of the complexity 
of tensor reduction for N > 5. However if one only needs to compute the 
rational part, it is not necessary to know the complete coefficients from tensor 
reductions. By the BDDK theorem [3], we know that many terms simply do 
not contribute to the rational part. Following this path of thought, the 
remaining problem is: is there an efficient way to compute these rational 
parts in one-loop QCD amplitudes directly from Feynman integrals? 

In this paper we show that there is actually a quite efficient and powerful 
method to compute the rational part directly from Feynman integrals. Be- 
cause we are concerned only with the rational part of the amplitude, there 
is no need for tensor reduction all the way down to scalar integrals. We only 
need tensor reduction to reduce the degree of the numerator by 2. So the 
original complexity of tensor reduction is bypassed in the computation of the 
rational part. 

In our approach of computing the rational part, we will exploit the theo- 
rem of Bern, Dunbar, Dixon and Kosower [3] and directly extract the rational 
part from the one-loop Feynman integrals. We also use the simple tensor re- 
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duction by using spinors as developed in [48, 66, 67]. We point out that 
the tensor reduction formulas used in our computations are actually quite 
simple, as one can see from eqs. (14) and (16) in Sect. 3. 

As we will demonstrate in this paper, the computation of the rational 
part is reduced to tree-level like calculations. As our method also applies to 
massive theory and theories with fermions, we envisage wider applications of 
our method in the computation of one-loop amplitudes, in combination with 
the D = 4 unitarity method [3, 48] and the efficient technique for computing 
generic unitarity cuts [38]. The once most difficult part of the one-loop 
amplitude can actually be attacked by the traditional technique. 

In this paper and the accompanying two papers [68, 69], we will develop 
our method and apply it to the computation of the rational parts of the 
one-loop 5- and 6-gluon QCD amplitudes. This paper mainly deals with the 
general theoretical formalism of the method. In [68], we show the efficiency 
of the method by computing the rational parts for the 5-gluon amplitudes 
for the two MHV helicity configurations by giving most of the intermediate 
steps. In [69] , we will present the results for the rational parts of the 6-gluon 
amplitudes for the two MHV and two NMHV helicity configurations. The 
rational parts of the 6-gluon amplitudes for the "split helicity" configurations 
are known already [34, 17]. Recently all one-loop maximally helicity violat- 
ing gluonic amplitudes were computed by Berger, Bern, Dixon, Forde and 
Kosower [70]. We refer the reader to [70, 69] for details about the explicit 
analytic results and comparisons. 

This paper is organized as follows: in Sect. 2 we set up our notation for 
spinor products and composite currents for sewing trees to the loop. Some 
simple tensor reduction formulas are given in Sect. 3. Starting from Sect. 4, 
we begin to develop the method of extracting the rational parts of Feynman 
integrals. We use the recursive approach to compute any Feynman integral 
as developed in [60]. In Sect. 6 and 7 we give explicit results for triangle 
and box integrals. In Sect. 8 we compute the correction terms to the naive 
D = 4 tensor reduction of box and triangle integrals, which arise from the 
ultra-violet divergent part of the box and triangle amplitudes. 
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2 Notation 

We mainly follow the notation of BDK [34] and the QCD-literature conven- 
tion for the square bracket [i j] . By abusing of notation the product between 
2 holomorphic spinors or 2 anti-holomorphic spinors is formed by a round 
bracket: 

(X i ,X j ) = (ij), (Aj, Xj) — (2) 

The scalar product between 2 vectors (written either in 4d vector notation 
or in 2 spinor notation) is also denoted by a round bracket. We have 

(X h X jl} X i2 X h ) = (X h , X i2 ) (X j2 , X h ) = (iii 2 )\j2ji\, (3) 
2 ki ■ kj = (Xi\, XjXj) = (ij) [j i}. (4) 

For spinor strings, we simply use (Xi\(k a + k b )\Xj) or (i\(k a + k b )\j] to denote 
(i-\(a + b)\j-): 

(XMK + kjfXj) = (i\(k a + k b )\j] = {i\(a + b)\j] 

= {i-\(a + b)\j-) = {ia)[aj] + {ib)[bj}. (5) 

We do not use gamma matrix traces. Instead we use bra and ket notation 
with multiple insertions of momenta: 

(i\hh--k n \j) = l <">[} " = odd . (6) 

j {i 1)[1 2] • • • {nj), n = even. 

Sometimes we also write (i\K\j) for with the understanding that 

sometimes the last j should actually stand for Xj and with the bracket ] . For 
example we have 

(ijA;! A; 2 Aj 3 |j] = <i 1> [1 2] <2 3> [3 j]. (7) 

For % — j the above notation is just the gamma matrix trace. For simplicity 
we will not write the slash: 

(i\ki k 2 k 3 \i] = ti-(ki k x k 2 k 3 ). (8) 

We note that the above notation only happens for an odd number of momenta 
inserted between 2 spinors (one holomorphic and one anti-holomorphic). Of 
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course the momentum can be either massless or a sum of several massless 
momenta. 

The sums of cyclicly consecutive external momenta are denoted gener- 
ically by Ki in a Feynman diagram. In our explicit computation we use 
k\2 = ki + &2 and ^234 = k 2 + k% + k^, etc. The kinematic variables are 
denoted as S12 = {k\ + k 2 ) 2 and s 12 3 = {k\ + k 2 + k^) 2 in a self explaining 
notation. For 6-gluon case we also have S123 = (£4 + k 5 + k 6 ) 2 by momentum 
conservation. 




(a) (b) (c) 

Figure 1: The composition of 3 external particles in tree amplitudes. 




(a) (b) (c) 




(d) (e) (f) 



Figure 2: The composition of 4 external particles in tree amplitudes. The 
blob denotes an expansion as given in Fig. 1. The explicit expression of 
Q(i+i)(j+2)(j+3) will not be given. 

For sewing trees to the loop, we define the following composite currents 
or polarization vectors: 

e i(j+l) = P\ e ii ki] ki + i) = — — — ■ M €j, ki + \) e i+ i 

- (ei+i, fci) Cj + ^(e;,e i+ i) (fcj - fci+i)), (9) 
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Q(i+l)(i+2) — P{£i(i+1), ki(i+l)] ^i+2) + -P(Q, e(i+l)(i+2), %+l)(i+2)) 

f (Q, £4+2) Q+l — ( e «> e «+l) e i+2 — ( e «+l? e i+2) €i) , (10) 

3 i(i+i)(i+2) V ^ 2 / 



+ 



where Sj(j + i)(j + 2) = (ki + /cj+i + k i+2 ) 2 . The above procedure is a simplified 
version of the general recursive calculation of the tree-level n-gluon ampli- 
tudes [71]. We note that e^+i) is anti-symmetric and e i ( i+1 )( i+2 ) is symmetric 
under the reversing of the order of the particles. This generalizes to com- 
posite currents with more legs, which we have not written down explicitly. 
The diagrammatic representations of Ci(i+i)(i+2) and ej(j + i)(j +2 )(i+3) are given 
in Figs. 1 and 2. 

In the formalism of [52] (see [53, 22] for reviews), the gluon polarization 
vectors are defined as 

where k is the momentum of the polarization null vector and q is the reference 
null vector. In terms of the holomorphic and anti-holomorphic spinors ( 
k = AA and q = rjfj), these polarization vectors can be recast as 



M)( t a) - ^° X f> e (-) (k . a ) - ^ Kf, P 

{r]X) , £ a ^k,q)- ^ 



In our notation, we will use e to denote the polarization vectors, and the 
relations between e and e is e = The troublesome y/2 will be absorbed 

in the overall coefficient of the amplitude. Given this notation, the bracket 
product of polarizations is more natural than the dot product : 

(ej(kj] qj), eT(k h g,)) = ej(kj] ■ e7(k h qi) = ^j^j ■ ( i3 ) 

In our calculations, it is also convenient to omit the denominators in the 
definition of the polarizations, and reinstate them in the last step. 



3 Tensor reduction of the one-loop amplitudes 

There is a vast literature on this subject. The original Passrino-Veltman 
approach [59] is quite general but it is not quite practical to obtain compact 
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analytic results. In fact the tensor reduction relations we will use for our 
calculations of the 5- and 6-gluon amplitudes are quite simple. It is based on 
the BDK trick [48] of multiplying and dividing by spinor square roots. To 
make more effective use of this trick, we have purposely chosen the reference 
momenta to make the tensor reduction simple. See [69] for further details 
about the specific choices of the reference momenta and the tensor reductions 
involved in the computation of (the rational parts of) the 6-gluon amplitudes. 



Figure 3: The tensor reduction for two adjacent same helicities generally 
gives three terms. 

There are basically only two different cases to consider. As shown in 
Fig. 3, the two polarization vectors have the same helicity. If we choose the 
same reference momentum (denoted by the spinor rj), we have 



propagators. The above tensor reduction formula is shown diagrammatically 
in Fig. 3, omitting the relevant factors. 

In deriving eq. (14), we assumed that p is a four-dimensional vector. 
Because pentagon and higher point one-loop amplitudes are ultra-violet con- 
vergent, the use of the above formula in tensor reduction is correct in dimen- 
sional regularization up to infinitesimal terms. However one must be quite 
careful to apply the above formula to the tensor reduction of the box and 
triangle tensor integrals because the difference is a finite rational part. Some 
correction terms must be included for tensor reduction with box and triangle 
tensor integrals. This also applies to the following tensor reduction formulas 
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where = (p + h) 2 , 1^ 



+ 
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given later in this section. We will compute these correction terms later in 
Sect. 8. 



2 



1" 





2 



2 



2 



Figure 4: For two adjacent same helicities, the tensor reduction for the com- 
bination of two diagrams is even simpler by a judicious choice of the reference 
momenta. 

An even simpler version of the above tensor reduction relation is to con- 
sider a combination of two diagrams together as shown in Fig. 4. The reduc- 
tion formula is: 



side of the above equations are read off directly from Feynman rules. These 
tensor reduction relations show quite clearly the simplicity of the diagrams 
when there are adjacent particles with the same helicity. 

For three adjacent particles with the same helicity gluons, we choose the 
following polarization vectors (omitting an overall factor for each polarization 





(p+fci)V(p-fc2) 2 ( P + hy( P - k 2 y 




vector): 



£4 — A5A4, 



£5 = 7?A 5 , 



^6 — A5A6. 



(17) 



Then we have 



e 45 = -7?A 4 + - 



1 (775) 



k 



(18) 



£56 = ??A 6 - - 



2 (45) 
1 fo5k 



(19) 



2 (6 5) 
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456 





(a) 



(b) 



(c) 



(d) 



Figure 5: For three adjacent same helicities, the tensor reduction for the 
combination of these four diagrams is also quite simple if we choose the 
reference momenta appropriately. 



and 

(e 45 ,p-k 45 )-±(e 4 ,e 5 ) = -(^P-M + ^H^-^), (20) 
(e 56 ,p-k 4 )-±(e 5 ,e 6 ) = faA 6 ,p - k 4 ) - \ M (J< B > - /«), (21) 

where = p 2 , 1^ = (p - A; 4 ) 2 , /( 6 ) = (p - A; 45 ) 2 and /« = (p - A; 456 ) 2 are 
inverse propagators. 

Now considering the three terms coming from the first 3 Feynman dia- 
grams shown in Fig. 5, we have: 

A 456 = (e 4 ,p- k 4 )(e 5 ,p- k 4 )(e 6 ,p- k 4 ) 

+ ((e45,P - ^45) - 2^ e 5 ))(e 6) P - ^4) ^ {5) 

+ (e 4 ,p- k 4 )((e 5(i ,p- k 4 ) - ^(e 5 ,e 6 ))/ (6) 

= ((A 5 A 4 ,p - k 4 )(r)\ 6 ,p - k 4 ) - (r)X 4 ,p- h)(\ 5 \ 6 ,p- k 4 )) I (6) 



2 (6 5) 



(A 5 A 4 ,p-A; 4 ) (7® -/W)/W, 



(22) 



by doing tensor reduction with /c 45 for the first term. This can be further 
simplified by expressing rj in terms of a linear combination of A 4 5. The final 
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result is: 

^456 



\w) {xiK p - h) J(4) J(5) + ^wi (As " 4 ' p - ki) /W /(6) 



1 {V 5) (A 5 A 6 , P - k 4 ) / W J( 5 ) + (A 5 A 4) p - fc 4 ) /« J( 6 ) 



2(4 5) 



[4 6] + 



2 (6 5) 
-(^Ai 5 i,p) 



(23) 



(4 5) (6 5) 

which has a nice symmetric property under the flipping operation 4^6. The 
last term in eq. (23) actually cancels the contribution from the last Feynman 
diagram in Fig. 5. 

For the case of different neighboring helicities, we can use the following 
reduction formulas: 

i \ ~ \ (\ \ \ (Xi\pk 2 Kk 1 p\fj) 
(Ai77,p)(AiA 2 ,p) = (2\K\1) 

= (2T^T) ( /W {Xl ~ X2,p) < 2 K^ + W 

+1^ (l\(p- h^K.k, - k 2 K,)\fj) - 1^ (l\pK 4 \l)[lf,] 
-/« (AiA 2 ,p) (2 1>[177] + (K 4 + h) 2 (AiA 2 ,p) (2 1> [1 77]), (24) 
(rj\pk 2 Kk l p\\2) 



(AiA 2 ,p) (t7A 2 ,p) = 
1 



(2\K\1) 
^Tfj(/ (1) (77|A; 2 ^3P|2> 

+/( 2 ) (771(^1 - k 2 K z ){p + fc0|2) - 1^ ( v \(k 2 + K 3 )|l)(AiA 2 ,p) 
+/ W (AiA 2 , p) (77 2) [2 1] - (k 2 + K 3 ) 2 (Ax A 2 , p) (77 2) [21]), (25) 

where J's are the various inverse propagators: 

/(i) = (p + fcO 2 , / (2) = p 2 , / (3) = (p - A; 2 ) 2 , (26) 



/( 4 ) = (p + fci + fT 4 



/(*) = ( p -k 2 -K 3 



(27) 



In eqs. (24) and (25), the momentum can be chosen as one of the nearby 
momenta to avoid the spurious pole associated with a composite momen- 
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turn. For two-mass-hard box, K can only be chosen as one of the composite 
momenta. 

Because of the complexity of the above tensor reduction formula for dif- 
ferent neighboring helicities, it is better not to use them directly. Luckily we 
are able to avoid using them directly for tensor reduction with 5- and 6-point 
diagrams by a judicious choice of reference momenta for all the polarization 
vectors. The details will be given in [69]. For two-mass-hard box integrals, 
we must use the above general tensor reduction in order to obtain compara- 
tively compact analytic expressions. We will use a slightly different reduction 
formula in Sect. 7.4 to compute the rational part of the two- mass-hard box 
integral. 

The above is just the first step for the tensor reduction. Of course this 
procedure can be applied recursively. By a rough examination of this recur- 
sive method, one immediately finds two problems: 1) the association of the 
resulting polarization vector with an external (composite) momentum may 
not satisfy the physical conditions; 2) the above formula is no-longer appli- 
cable if one has a massive external momentum (a composite one arising from 
the reduction or a pinched line by sewing the tree to the loop). There are 
other methods to do further tensor reduction, some involving Gram deter- 
minants. Exactly because of these problems, tensor reduction is usually the 
most difficult part and the bottleneck for directly computing the one loop 
amplitude. Quite elaborate methods are developed to tackle these problems. 
See, for example, [8, 61]. 

Because of the complexity of doing further tensor reduction, we immedi- 
ately see the problem why directly computing the amplitude by using Feyn- 
man rules is an extraordinarily difficult task for higher point amplitudes. 
There are mainly two difficulties to overcome: too many diagrams and the 
complexity for tensor reduction (especially at a later stage of tensor reduc- 
tion as mentioned in the above). By using computers it is not too difficult 
to manage the numerous Feynman diagrams (of the order 1000). But the 
complexity of tensor reduction with the appearance of spurious poles is ac- 
tually the bottleneck for analytic computation. Tensor reduction is also the 
bottleneck for doing numerical calculations. See, for example, [62]. 

In contrast we also see that why it is possible to compute the rational 
part by using the conventional Feynman integrals. First by using the super- 
symmetric decomposition, the number of Feynman diagrams is about 50 (1 
hexagon, 6 pentagons, 15 boxes, 20 triangles and 15 bubbles) by only com- 
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puting the scalar loop contributions. Second by computing only the rational 
part, it is not necessary to do tensor reduction all the way down to scalar 
integrals. One needs only to do tensor reduction to reduce the degree by 
2 (see next section). So tensor reduction does not complicate the analytic 
expressions significantly. In fact for special helicity assignments, there is an 
almost mutual cancellation between higher point diagrams and lower point 
diagrams, as we demonstrated in eqs. (14) and (16). This is a manifestation 
of gauge invariance. This property can be used to check and to organize the 
results of our calculation. It is very important to have some "local" can- 
cellations before adding all the results together in order to obtain relatively 
compact analytic results for the rational parts of QCD amplitudes. 

4 The BDDK theorem and the structure of 
one-loop amplitudes 




Figure 6: A generic one-loop diagram with external momenta K±, • • •, K m . 
p is the internal momentum between the external lines K\ and K m . 

A generic m-point one-loop Feynman diagram shown in Fig. 6 and its 
integral is given as follows (by using Feynman rules): 

OfiP)} = J —if 2 p2{p _ Ki y P \ {p ~ ^ (28) 

where f(p) is a polynomial function of the internal momentum p. For phe- 
nomenologically interesting models and by choosing a suitable gauge, the 
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degree of f(p) is always not greater than m. f(p) also depends on the exter- 
nal momenta k (K^s are sums of cyclicly consecutive external momenta fc's) 
and the polarization vectors tj. For f(p) = 1 it is called the scalar integral. 
The strategy of computing I^[f(p)] is to reduce it recursively into lower de- 
gree polynomials and/or lower point integrals. It is a well-known result that 
In If] can ^ e generically written as: 



up to infinitesimal terms arising from tensor reduction from pentagon or 
higher point diagrams. Here rational functions of the external mo- 

menta when all polarization vectors are written in terms of spinor products. 
We note that these coefficients also depend on the (arbitrary) space-time 
dimension D in dimensional regularization (in the FDH scheme [72]). 

A brute-force computation of these coefficients from Feynman integrals 
is an impossible task for 6 or higher point amplitudes. The 5-point case 
was computed by string-inspired method by using a table for all Feynman 
parameter integrals (see below) [7]. However the string-inspired method is 
still not powerful enough to compute even the 6-gluon amplitude due to the 
complexity of the Feynman integrals and the intermediate expressions. 

For physical application what we actually need is an expansion in e of 
the above formula for D = 4 — 2e, up to finite terms. If we forget the 
infrared divergence for the moment, there are only simple pole (-) terms in 
the scalar integrals Jfg 2 [1] with rational coefficients. So we can write the 
n-gluon amplitude (addition of all contributing I^[f(p)] from all Feynman 
diagrams) as follows: 



i i 

+ ]T c 2 ,i(e, k; 4) J 2 DW [1] + (rational function) + 0(e). (30) 



For supersymmetric theories, Bern, Dixon, Dunbar and Kosower [3] proved 
a theorem which states that the rational function is exactly zero. This is due 
to the better ultra-violet behaviour of one-loop amplitudes in supersymmet- 
ric theories. So we need only to compute the rational coefficients (called the 



Of] 



+ £ c 3 ,,(e, k; D) + £ c 2 ,(e, k; D) 



(29) 
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cut-constructible part hereafter) exactly at D — 4. What they proved is ac- 
tually a more general theorem: if f(p) is a polynomial (in p) of degree m — 2 
or less, the rational part for I®[f(p)} arising by expanding in e is exactly 
zero. Generally speaking, the rational part is non-vanishing for degree m 
and m — 1 polynomials. 

For a non-supersymmetric theory like QCD we also need to compute the 
rational function (called the rational part hereafter). In a series of papers 
[3, 48, 50], Bern, Dunbar, Dixon and Kosower have developed a method of 
computing the cut-constructible part, i.e., the coefficients Cj^(e, k; 4), from 
4-dimensional unitarity. The nicety of 4-dimensional unitarity is that all 
the ingredients in the unitarity relation are on-shell quantities. However 4- 
dimensional unitarity loses all information about the rational function part 
in eq. (30). One must use some other methods to compute the rational 
part. In [48], they use the factorization properties with trial and error. The 
difficulty of computing the rational function prevents the wider application of 
the unitarity method to calculate more general amplitudes. As we mentioned 
in the introduction, the rational function part could be computed by going 
to /^-dimensional unitarity [50, 42]. 

In our papers, the precise definition of the rational part R n of the QCD 
amplitude is as follows. The n-point one loop color-ordered partial amplitude 
A^i defined in eq. (2.10) in the second paper of [34] can be decomposed as 

A nl = ^(E^(e^;4)/ 4 DW [l]+Ec3,(6^;4)7 3 I3W [l] 

i i 

+ E^(6,A;;4)/ 2 DW [l]+2i? n + 0( e )), (31) 

i 

where R n is the rational part for one real scalar circulating in the loop. We 

will use notations like R( h + ++) to denote R n in different cases in [68] 

and [69]. 

In the following sections we will exploit the BDDK theorem to compute 
the rational part directly from the Feynman integrals. By using the recursive 
relations satisfied by the tensor integrals we will derive the recursive relations 
for the rational parts by making an expansion in e. Our integration method 
of computing the rational part may also be used to compute the rational part 
by using D-dimensional unitarity. We use then recursive relations to derive 
explicit formulas for the rational parts of all bubble and triangle integrals 
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in Sect. 6. In Sect. 7, we derive the formulas for box integrals up to two- 
mass-hard boxes. Formulas for 3-mass and 4-mass box integrals can also be 
derived. They are much more complicated than the two-mass-box formulas. 
Fortunately they are not needed in the computation of 6-gluon amplitudes 
and will not be given here. We note that the recursive relations for tensor 
integrals and the rational parts can also be derived for massive internal loop 
and/or external fermion lines. For simplicity all formulas are given only for 
the cases with vanishing internal masses. 



5 The recursive relations of one-loop ampli- 
tudes 

In this section we study the recursive relations of one-loop amplitudes [60, 
61, 8]. By using Feynman parametrization we have 



d u p 1 

i n D/2 p2^p _ fc^y . . . (jj _|_ 

1 )"!•(// D/2) [d n a 5(1 " >>;i 



(a-S-a) n -~ 



where 



and the matrix S, 



a • S • a = ^ cij Oj Sij 



( k\ {h + k 2 f 
* 



s = 



k\ 



* 
* 



(h + k 2 + ■ ■ ■ k n ^) 2 \ 

(k 2 + k 3 + ■ ■ ■ k n ^f 





* 







/ 



(32) 



(33) 



is an n x n symmetric matrix of external kinematic variables (extension to 
massive loop is straightforward). For tensor integral In[f(p)} it is given 
by the Feynman parameter integral with an extra polynomial of a in the 
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numerator 2 : 



I°[g(a)} = T(n-D/2) d n a 



6(1 - E< a«)g(a) 
(a-S-af^ 



(34) 



The degree of g(a) is the same as the degree of f(p) in p. 

As explained in [60] , a set of recursive relations for these tensor integrals 
can be derived by performing the following integration: 



F = T(a) 



1 rl— a-n-l r\—a\—a,2- 

da„_i / da n _ 2 ■ • • / 

o Jo Jo 



-a„_i 



da r 



x 



9 



<9a r , 



f(a) 



(a-S- a) a 



= 1 — 0.1 On- 



(35) 



in two ways: one by partial integration and one by direct differentiation. 
Then we obtain the following recursive relation: 

-2r(a + l) fd^l-aj^M^M 
v 'Jo (a • S ■ a) a+l 



Jo (a • o ■ a 

-F(a) C d n a8(l-a) — 
Jo (a 

-r(a) ^(TaSn-a) 
Jo 



a„=0 



S ■ a) 
d m f(a) - d n f(a) 



(a- S ■ a) c 



(36) 



By carefully examining the composition of a • S • a for a rn = one recognizes 
that this corresponds to a pinching limit of the original Feynman diagram: 
k m -i, k m — > k m _i + k m . By setting a = n — 1 — y and using the definition 
for I^if], the above recursive relation translates into the following form: 



I°l-2f(a)((S-a) m -(S-a) n )} 



&\f(a)\ ~ &\f(a)\ - In +2 [d m f(a)] + I° +2 [d n f(a)]. (37) 



p£»(m) 



2 We note that the difference between Feynman parameter integral I n [g(a)] and / n [/(p)] 
is only the (— 1)" factor. For even n and f(p) — g(a) — 1, they are the same and we will 
not distinguish them in this case. 
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Here In-i[f( a )] denotes the [n— l)-point one-loop Feynman integral obtained 
by pinching and ki. Coupled with one more equation from the delta 
function: 

n 

E4 D [/(aH = /£[/(«)], (38) 

we can solve if?[f(a)ai] in terms of in-T\f(o)], ^n +2 [9 m f{a)} and l£[f(a)]: 
In[f(a)a t } = ^^[/(^l+^/^^/C^l+^^f/Ca)]. (39) 

The above reasoning goes through so long as we can invert the relevant 
matrix. This matrix is 



—2(Sij — S n j) 
1 ••• 1 



(40) 



and it is singular for n > 6. For our purpose of computing the rational part 
we will only need these recursive relations for n < 4. Higher point tensor 
integrals are reduced directly in D = 4 as we have done in the last section. 
For further discussions about the tensor reduction and its close relation with 
the above recursive relation, we refer the reader to [61, 8]. 

The above recursive relation is not symmetric for all a^. The symmetric 
recursive relation can be obtained by firstly solving I®[f(a)(S ■ a) n ] and 
substituting it back to the system of equations. By taking / = gi(a) a m (gi(a) 
is a homogeneous polynomial of degree / in a), we can solve I^[gi(a) (S ■ a) n ] 
to get: 

2fffo(a) (S ■ a) n ] = in-i[gi{a)\ + {n-l-l-D) i° +2 [ gi {a)\ + i° +2 [d n9l {a)\. 

(41) 

Substituting this back into eq. (37) and multiplying both sides of the equation 
by and then summing over m, we have [61]: 

inh(a)a t } = ±(n-l-l-D) lt i^[ gi (a)] 

+ ^E^^t^Wl + ^E^ 1 ^ 2 ^^)], (42) 



where 



7i = E^- ( 43 ) 
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One can check that the above recursive relation is equivalent to the following 
recursive relation: 



#W(a)] = Pa (4 D(j) [/(«)] + 4 D+ W(a)]) + | ff[/(a)], (44) 

Pa = \( s *- 1J g)> A = ^7, (45) 

The good point of this recursive relation is that all coefficients have no explicit 
dependence on the space-time dimension D and so it is well suited to compute 
the rational part. 

We note that the above reasoning would go through if S is invertible. In 
the two-mass and one-mass triangle cases, S is not invertible. Nevertheless 
we can still use the above formulas by taking a limit from the general 3-mass 
triangle case. The point is that all we need is to have a well-behaved limit for 
the quantity and ^. We have verified that directly taking the massless 
limit gives the same result as one would obtained by solving the non-singular 
system of equations (37) and (38). 

From eq. (44) we see that higher dimensional integrals also appear in the 
recursive relations. These higher dimensional tensor integrals can be reduced 
to even higher and/or lower point tensor integrals by using these recursive 
relations repeatedly. At the end only scalar integrals are left. These higher 
dimensional scalar integrals can be reduced to lower dimensional and lower 
point scalar integrals by using the following recursive relation: 



n-l -D) A 



(46) 



This equation is derived from eq. (42) by setting gi(a) = 1 and summing over 
i. The explicit dependence on the space-time dimension D in eq. (46) is very 
important. Otherwise all the rational coefficients in eq. (29) would have no 
explicit dependence on D and there would be no rational part. 

To be more specific, we give some explicit examples for bubble and tri- 
angle integrals in what follows. 

For bubble integral we have 
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where K is the sum of momenta on one side of the bubble diagram. For 
K 2 = this integral is in dimensional regularization. So we will assume 
K 2 7^ hereafter. By direct computation we have the following results 
(D = 4-2e): 



e(l -2e) 



2\-e 



r(i + e)r 2 (i- e ) 

r(l-2e) 



/f[a 2 ] = J 2 > 2 ] 



2-e 



2(3 - 2e) 

The recursive relation (46) becomes: 



if 2 

2( J D-1) /2D[1] - 



(48) 

(49) 
(50) 



(51) 



This can be applied recursively to compute arbitrarily higher dimensional 
bubble integrals. 



P-K, 




P + K 3 



Figure 7: A generic three-mass triangle diagram. The 2-mass triangle dia- 
gram is obtained by setting one of momentum to be massless, for example 
K x = ki and k 2 = 0. 

A generic triangle diagram is shown in Fig. 7 and the integral is 3 

I?[f(p)\ = J p2 {p _J^\ p + K3) 2- (52) 

3 A minus sign is not included in the definition of I^[.f{p)\- 
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This integral is finite (free of ultra-violet (for degree 3 or less polynomial 
f(p)) and infrared divergences) for generic external momenta (i.e., Kf ^ 0, 
% = 1,2,3). The explicit formula can be found in [3]. For degenerate cases 
we have 



t d U] i r(i-e) 2 r(i + 6) (-if 2 2 )--(-iff)- 2 

h [1] = ~7 r(i- 2e ) (-kI)-(-k!) ' ^ = °' (53) 

= - ^^ra-L^ 6 ^ -^ 2 )" 1 ^ *? = *2 = o. (54) 

In the next section we will use these results and the more general recursive 
relations to derive the rational parts of the bubble, triangle and box integrals 
by making an expansion in the space-time parameter e. 

6 The rational parts of the triangle integrals 

6.1 e-expansion and the rational parts: the bubble in- 
tegrals 

The e-expansion of the scalar bubble integral is: 

I 2 D [1] = - e +0(l). (55) 
By using this result in eq. (51) we have (D = 4 — 2e): 

/f +2 [l] = ^/ 2 D [l]+^ + 0(e), (56) 

/ 2 D+4 [i] = {J ^ i ^+ 4 -W- +o{e) - (57) 

We note that the first term on the right-hand of the above two equations still 
depends on D through the scalar integral /^[l]. What is important is the 
second term which is a pure rational function and does not depend on the 
space-time dimension D. This is the rational function part we need to keep 
track of later. 

By making an expansion in e for eq. (50), we have: 

i 2 D [al} = I 2 D [at} = \ I 2 D [1] + 1 + 0(e). (58) 
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From this equation we read off the rational part of i^t ?] and i^a^] as ^. 
Translating back to the momentum integral the result is: 



+ ^((fi.A')(e2.A')-2A' 2 (6 1 ,e 2 )). (59) 

This result has already appeared in [50] (eq. (31) on p. 133). We interpret 
the second term in the above equation as the rational part. If we discard the 
first term we can simply write: 

I 2 [(e u p) (e 2 ,p)] = 1 (( Cl , K) (e 2 , K) - 2 K 2 ( Cl , e 2 )), (60) 

by dropping also the explicit dependence of I 2 on D. However we still retain 
this dependence for higher dimensional Feynman integrals and simply drop 
all the cut-constructible parts. Explicitly we have: 

i 2 [l] = 0, (61) 
i 2 [a\} = i 2 [al] = -i 2 [ ai a 2 ] = ^ (62) 

I™[1] = ^. (63) 

6.2 The rational parts of the higher-dimensional scalar 
integrals 

For higher dimensional scalar integrals we can use the recursive relation 
eq. (46) to derive their rational parts. 

For three-mass triangle the explicit recursive relation is: 

?3m(D+2) m 1 [ f3 m . , . x r (l) 

h [1] = (2 - £)) A [ 2siS2S3/ 3 + s 2{si + s 3 - s 2 ) P 2 ' 

+ s 3 (si + s 2 - s 3 ) I { 2 2) + si (s 2 + s 3 - si) 4 3) ] , (64) 
A = s 2 + §2 + S3 - 2(si s 2 + si s 3 + s 2 s 3 ), (65) 

where Sj = iff, i = 1, 2, 3. 
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For 3- mass triangle integral 1% [1] is regular as D — > 4. By using the 
e-expansion of we have: 



^(D +2)[l] = __L [2 aia2S 3/*» + S2 ( ai + a3 - a2 )/l 



3 3m + S 2 (Sl + S3-^ W{1: 



+s 3 (si + s 2 - s 3 ) # ; + si (s 2 + s 3 - si) J^J + l - + 0(e). (66) 

So the rational part of (the Feynman parameter integral) J 3 D+2 [1] is |. This 
result also applies to the two-mass and one-mass triangle, although there are 
intricacies of infrared divergences. This can be explicitly checked by using the 
following explicit recursive relations for the two-mass and one-mass triangle 
integrals: 

/3 m(D+2) [l] = ( 4 3) - IP) , (67) 

3 1 J D-2 \ Sl -s 2 2 s 1 -s 2 J y ! 

H m{D+2 \l] = ^ 4 3) , (68) 

and the explicit formulas for J^^'s. We note that the above recursive re- 
lations can be derived either by taking the limit s± = K\ — > (2-mass) 
or Si 5 2 = 7Ti2 ~~ *■ (1-mass) or by solving eqs. (37) and (38). In above 
eqs. (67) and (68), the terms with infrared divergences really do not appear. 
This is because the coefficients of the infrared divergent terms in eq. (65) 
are automatically zero. The infrared divergences also do not contribute to 
the rational parts in the four point cases, which can be checked by explicit 
calculations. A heuristic argument is that since higher dimensional integrals 
are free of infrared divergences, we expect that there should not be infrared 
divergence in the square bracket on the right hand side of eq. (46) like in eq. 
(67) and (68), and hence the infrared divergences do not contribute to the 
rational part. 

The end result of the above analysis is 

J 3 D+2 [1] = \ + c(s) + 5>(*)/? (<) [l] + 0(e), (69) 
z i 

where c(s) and q(s)'s are (rational) functions of the external kinematic vari- 
ables Si and the polarization vectors. This formula applies to all possible 
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triangle integrals. This shows explicitly that the rational part of the higher 
dimensional scalar integral J.f +2 [1] is a purely ultra-violet effect. 

Setting the cut-constructible part to 0, we can effectively write the fol- 
lowing formula for the rational part: 

4 D+2 [1] = \ (70) 

We have also studied in detail the box integrals by using the explicit 
formulas of [3]. We checked all the degenerate cases. The results can be 
simply stated as follows: 

Jf+ 2 [1] = /f +2 H=0, (71) 

4 D+4 [1] = 4 (72) 

by dropping all the cut-constructible part. For 5-point integrals we have 

/f+2[1] = /^ [ai]= /^ [aia .] =0 , (73) 
f D +4[1] = j 5 ^4 W=0; (?4) 

*™ - ^ (75) 

although they are not used in the computation of the QCD amplitudes. 

The rational part for the tensor integral is computed by first transforming 
it into Feynman parameter integrals and then using the recursive relations. 
The recursive relation for the rational part is exactly the same as for the 
complete Feynman integral. However all lower degree (m — 2 or less for m- 
point) Feynman integrals can be set to zero by BDDK theorem. Effectively 
the recursive relations are truncated. In the next few sections we compute 
explicitly the rational parts for triangle and box integrals. From hereafter 
all Feynman integrals will mean their rational parts, except explicitly stated 
otherwise. 



6.3 The triangle integrals: the general case 

For triangle integrals we will consider 2 cases: the two-mass triangle and 
the three-mass triangle, although the two-mass case can be obtained simply 
by setting one of the masses to be 0. The reason for doing this is that 
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the formulas simplify greatly for the two-mass triangle integral. For some 
computations one may only need the simplified formulas (for example the 
5-gluon amplitudes in [68]). We will first give the formulas for the general 
three-mass triangle integrals. 

The three-mass triangle diagram is shown in Fig. 7 and the external 
momenta are denoted by Ki (i = 1, 2, 3). We use the convention of denoting 
a light-like momentum by a lower case k, i.e. k 2 = 0. 

We will give the rational parts for both the degree 3 and degree 2 polynomials 4 . 
First we make the following definitions: 

w c x _ f ifiuP) (e 2 ,p-K 1 )(e 3 ,p + K 3 ) 

/3(ei ' e2,e3) = pUp-Ktffr + Ktf ' (76) 

T ( e f \ = f ( £ ^P) (77) 

31 u 2) ~ J p 2( p _ Kl y( p + K3 y V) 

The actual computation of the rational part is done by using Feynman 
parametrization and we have: 

h(ei,e 2 ,e 3 ) = -I 3 [(e 1 ,a)(e 2 ,a)(e 3 ,a)] 

+ f>,e m )/ 3 D+2 [( ej+2 ,a)], (78) 
i=i 

h(ei,e 2 ) = -I 3 [(e 1 ,a)(e 2 ,a)] + ^(e 1 ,e 2 ), (79) 

by using the previous result for ig 2 [1]. In the above we have used the 
following shortened notation: 

(ei,a) = (e!, Ki) a 2 - (e!, K 3 ) a 3 , (80) 
(e 2 ,a) = {e 2 ,K 2 ) a 3 - (e 2 ,K 1 ) a u (81) 
(e 3 ,a) = (e 3 , K 3 ) a x - (e 3 , K 2 ) a 2 , (82) 

where Feynman parameters as used in Sect. 5. In order to give 

compact formulas for the various quantities appearing in eqs. (78) and (79), 
we first define the following functions: 

-^0(51,^2,53) = --^10203] 



4 It is also possible to derive the rational parts for higher degree polynomials but they 
have no practical usage in application to computations in the electroweak theory and QCD. 
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10sis 2 s 3 (si + S 2 + S3) 
3A 2 6A 
^i(si,S2,s 3 ) = -h[a 1 a 2 (a 2 + a 3 )] 



^2(51,52,53) 
^3(51,52,53) 



3A 2 


+ 


3A 


-/ 3 [a 2 a3(a3 + ai)] 






5(5 2 + s 3 -si)s 3 Si 


+ 


(5 2 - Si) 


3A 2 


3A 


-/ 3 [a 3 ai(ai + a 2 )] 






5 (s 3 + si - s 2 ) sis 2 


+ 


(S3 - s 2 ) 



3A 2 3A 



where Sj = K 2 and A = s 2 + s 2 + s| — 2(sis 2 + s 2 S3 + s 3 si) is the 
determinant" for the triangle diagram. 
By using these functions we have 

J 3 [(ei, a) (e 2 , a) (e 3 , a)] = F (si, s 2 , s 3 )((ei, Xi) (e 2 , Ki) (e 3 , X 2 ) 
+( Cl , AT 3 ) (e 2 , tf 2 ) (e 3 , tf 2 ) + (e,, K 3 ) (e 2 , ^) (e 3 , tf 3 ) 
+(ei, AT 3 ) (e 2 , Xi) (e 3 , K 2 ) - (ei, #1) (e 2 , tf 2 ) (e 3 , tf 3 )) 



+ XX e i> ( 6 2> ( e 3, Ki) Fi(si, S 2 , S3) 

1=1 

+i((5i - s 2 - s 3 ) (ei, ifi) (e 2 ,K 1 ) (e 3 ,K 3 ) 



2A 

+ (s 2 - s 3 - si) (ei, Ki) (e 2 , X 2 ) (e 3 , X 2 ) 
+ (S3 - si - s 2 ) (ei,X 3 ) (e 2 ,K 2 ) (e 3 ,X 3 )), 



lf«[( ei ,a)| = fo.jr.-jr,) (1 + 52(83 

+(elIA . J) ( flzf ^+ flzf! ) I 

+K A3) . 
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h D+2 [(e 3 ,a)} 




) 



(90) 



and 



= ± ( Sl (( Ci , tf 2 ) ( Cj , tf 3 ) + (q, tf 3 ) (ej, ^2)) 
+ s 2 (( £i , tf 3 ) (e„ + (ei, Xi) (e„ X 3 )) 
+s 3 ((ei, Xi) (e„ tf 2 ) + ( Ci , tf 2 ) (e„ XO)) . 



(91) 



We note that some formulas (and also the function Fi) in the above are 
related by permutations. We purposely wrote down the complete formulas 
in order to see the pattern. 

6.4 The triangle integrals: the two-mass case 

For two- mass triangle we set K\ — k\. Then we can derive the simplified for- 
mulas in the two-mass triangle case by setting s\ = in the above formulas. 
Also we redefine (e 3 , a) to be (e 3 , k\) a 2 — (e 3 , K 3 ) a 3 . This corresponds to the 
change of (e 3 ,p + K 3 ) to (e 3 ,p). This gives a more symmetric form for the 
result, as one can see from the Feynman integral representation by doing the 
transformation 2^3. We list all the explicit formulas here because they are 
used heavily in the actual computation of the 6 particle amplitudes. This 
is necessary to obtain compact analytic formulas for the rational part of the 
amplitude. We have: 



h[(ei,a) (e 2 ,a) (e 3 , a)] 




7 v. ( 6l > K ^ ( e2 ' ( e3 ' 

(«2 - S3) 2 

((e 2 ,k 1 )(e 3 ,K 3 )-(e 2 ,K 2 )(e 3 ,k l )). (92) 



6(s 2 - S3) 




(93) 



; L>+2 
3 



[(e2,a)] 



3g (c2,A;i) + - (e 2 ,X 2 ) 
5 (es,fci) -Ues,K 3 ) - 



(g2 + S3) (f2,fcl) 

12(s 2 -s 3 ) 



(94) 



4 D+2 [fe,a)] 




(95) 



12(s 2 - s 3 ) 
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7 3 [(e i ,a)(e j ,a)] = - ^ * ^ (e j; fci) (e,, fci) 

(fa, fci) (e 3 -, Ji~ 2 ~ gg) + fc i) (fr, ~ ^3)) , g6 ^ 
4 (s 2 - s 3 ) 

We note that in the last formula the double pole term is absent if one of the 
polarization vectors is associated with the first momentum k\ and satisfies 
the physical condition (e, fci) = 0. 

By using the above results we have 



^1^2) = / 7 
J 1 



d D p (e 1 ,p)(e 2 ,p) 



i7T D/2 ptfy-ktfip + Ks) 2 

1 . . {Kl + Kl) . , w 

= 2 (Cl ' Ca) + 2(4-4) a (£l ' fcl)(£2 ' fcl) 
((ei, K 2 ) (e 2 , fci) - (ei, fci) (e 2 , Jig)) 

2(Ki-Kl) ' ^ 



1 , x . ('I - /\' 2 ) ('2- /•'! ; , , „ 

Wei, ^2) = - (ei, e 2 ) + yjqZTJq) ' ^ = °' ; 



and 



d D p (ei,p) (e 2 ,p- fci) (e 3 ,p) 



iTr^/ 2 p 2 (p - fcO 2 (p + X 3 ) 2 
^ ((e 2 , 4 X 2 - 7 fci) ( ei , e 3 ) - (2 <- 3) + 4(e 1? X 2 ) (e 2 , e 3 )) 
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(€!,K 2 ) (e 2 ,fci) (e 3 ,fci) 



6(K|-^ 3 2 ) 5 
(ei, tf 2 ) ((e 2 , fci) (e 3 , # 3 ) - (e 2 , #2) (e 3 , fci)) 



+ gf) 
12 (X| - Xf) 



6 (Xf - Xf) 

((e 1 ,e 2 )(e 3 ,fci) + (e 1 ,e 3 )(e 2 ,fc 1 )). (99) 



The above formula is anti-symmetric under the exchange 2 <-> 3 by not- 
ing (ei,/\~ 2 ) = — (£1,-^3)- For later application we denote the two-mass ra- 
tional part of 73 by w h ere { denote the massless external line. To 
distinguish the two possible two-mass triangle diagrams with the same mass- 
less external line in the 6-gluon amplitude case, we put a tilde on if" 1 for 
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one of the two-mass triangle integrals with 1, 3 and 2 external momenta 
in a clockwise direction. Referring to Fig. 7, J 2m ^ has external momenta 
{ki,k i+ i + k i+2 ,k i+ 3 + k i+4: + k i+5 }, whereas J 2m ^ has external momenta 

7 The rational parts of the box integrals 

A generic box diagram is shown in Fig. 8. The kinematic variables s and t 
are defined by following the standard notation: 

s = (K, + K 2 f = (K 3 + K 4 )\ t = (K 2 + K 3 ) 2 = (K A + K x f . (100) 

For our purpose of computing up to 6-gluon amplitudes, we need to consider 
up to two-mass boxes. There are two kinds of two-mass boxes: the two- 
mass-hard box and two-mass-easy box. By a judicious choice of reference 
momenta the two-mass-hard box does not show up in the computation of 
MHV amplitudes. So we will discuss only the two-mass case here and set 
K\ — ki, etc. by following the convention of writing the light-like momenta 
in lower case k. The one-mass box case is obtained as a special case of the 
two-mass-easy box case by setting further K 2 = k 2 . We will list the explicit 
formula for one-mass box for quick reference and the ease of use. 




Figure 8: A generic box diagram, p is the internal momenta between K 4 and 
K\. Other internal momenta are also shown explicitly. 
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7.1 The box integrals of degree 3 polynomials: the 
two-mass-easy case 

Generally we need to compute the rational part of the following box integral: 



(e2,p) (e 3 ,p) 



i7T D/2 p 2( p _ Ki y^ p _ Kl _ K 2 y( p + K4 y 



(101) 



The complete rational part in the general case is quite complicated and 
should be avoided. From our experience it is always the case that at least 
one of the polarization vectors satisfies the physical condition for one of 
the massless external momenta k\ or k 3 . In fact one can always expand 
an arbitrary 4-dimensional vector in terms of the 2 independent spinors of 
the two external massless momenta. So we can assume that ei satisfies the 
physical condition: (ei,ki) = 0. To be specific we can take ei = 77 Ai. The 
negative helicity case can be obtained from this case (the positive helicity one) 
by conjugation. If one of the polarizations satisfies the physical condition for 
&3 we can rotate the two- mass-easy box diagram by tt and relabel k 3 as k\. 

By explicit computation, we found that if the reference momentum of ei 
is k 3 , the rational part becomes quite simple and is given as follows: 



/4(A 3 Ai, e 2 , 63) 



<3|tf 2 |l] 



(£2,^3) (^3,^3 



{t2,h) (e 3 , h 



{Kl-t){Kl-s) (Kl-s)(K!-t) 



(102) 



By using this result, the computation of the rational part of the degree 3 
polynomial can be proceeded by changing the reference momentum of e\ to 
k 3 . This is equivalent to expanding the spinor in terms of A 1,3: 



, = M Al + M A3 . 



(13) 



(31) 



(103) 



In so doing we also generate 2 triangle diagrams which have been computed 
in last subsection. Explicitly we have: 



^4(^1,62,63) = 



+ 



(?7 1) (3|if 2 |l] 
(31) 2 

fa 3 ) r?2m, 



(€2, ^3) (^3, ^3 

[Kl-t)(Kl- 



(€2,fcl) Ml) 

(Kl-s)(Kj-t) 



(13) 



£2,63) 



/ 3 2m (e2,e 3 ) 



(104) 
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I~i m (e 2 ,e 3 ) = ]- (e 2 , 6 3 ) + igf + 1 (e 2 ,fc 3 ) (e 3 , fc 3 ) 



+ 



2^'^ 1 2(K 2 2 -t)2 

((€2, fcg) (£3, gg) ~ (^3, fc 3 ) (£2, ^4 + fci)) 

2(KI - 1) 



l! m (^,es) = \ (e 2 , 63) + 2 ( ( ^ ± ^ (e 2 ,fc 3 ) (e 3 ,fc 3 ) 

((e 2 , fc 3 ) (e 3 , if 4 ) - (e 3 , fc 3 ) (e 2 , fci + K 2 )) , . 

+ 2(Xf - S ) • (1 ° 5) 

Taking into account the anti-symmetric property of the product (A3A1, i^ 2 ) = 
— (A3A1, K4), the above formulas are actually symmetric (which must be the 
case) under the interchange 1 <-> 3 and 2^4 (the polarization vectors 
are kept fixed because they should be invariant under the permutation by 
themselves) . 

7.2 The box integrals of degree 3 polynomials: the 
two-mass-hard case 

For the two-mass-hard box case, we follow the same strategy. For e\ = AiA 2 
or A 2 Ai, we have 

J 4 2 ^(A 1 A 2 ,6 2 ,6 3 ) = { -^f^h(e 2 ,e 3 ), (106) 
/4 2m/l (A 2 A 1)e2 , e3 ) = ^H/ 4 ( e2 , e3 ), (107) 

where 

h(e 2 ,e 3 ) = (e 2 ,h)(e 3 ,K 4 ) + (e 2 , if 4 )(e 3 , fci) + (e 2 , k 2 )(e 3 , K 3 ) 

+ (e 2 , K 3 )(e 3 , fc 2 ) - 1 [2 {KIkI - 1 2 + 5)(e 2 , fc 12 )(e 3 , fc i2 ) 

+ (X 3 2 + K 4 2 - s - 2 ((^ 3 2 - K\ + s)(e 2 , K 4 )(e 3 , K 4 ) 
+(K 2 4 -K 2 3 +s) (e 2 ,K 3 )(e 3 ,K 3 ))~ 

+ ^^(e2,fci)(e 3 ,fci) + ^^(e 2 ,fc 2 )(e 3 ,fc 2 ), (108) 
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and 

5 = K*Kl-(K* + Kl)t+(s + t)t, (109) 

A = A(k 2 12 ,KiKl), (110) 

A(si,s 2 ,s 3 ) = si + S 2 + 4 - 2(s 1 s 2 + s 2 s 3 + S 3 Si), (111) 

are functions of the external momentum invariants. In particular A is the 
Gram determinant of the three-mass triangle integral arising from the tensor 
reduction of the two-mass-hard box integral. 

The general case is obtained by changing the reference momentum: 

(77 1) 

^4(^1,62,63) = — — / 4 (A 2 Ai, e 2 ,e 3 ) 



+ 



(21) 

M 

(12) 



Ii m (e 2 ,e 3 )-Ii m (e 2 ,e 3 ) , (112) 



1/ \ (Kl + t) 



irises) = - (62, e 3 ) + 2 V I _ I (6 2 , h) (63, k 2 ) 

((€2, fc 2 ) (63, if 3) ~ (fa, fa) (eg, -K4 H- fci)) 

and 7| m (e2,63) is the three-mass triangle integral with external momenta 
{ku, K 3 , K4}. There are two triangle diagrams from this reduction and one 
is a three-mass triangle diagram. 

7.3 The box integrals of degree 4 polynomials: the 
two-mass-easy case 

Now we discuss the computation of the rational parts for degree 4 polyno- 
mials. First we define 

T , f d D p (ei, p) (6 2 , p- K l )(e 3 ,p- K l2 ) (e^p + Kj) 

4l61 ' C2 ' C3 ' 64J " J m°l 2 p 2 (p - Xi) 2 (p - K 12 f{p + K 4 y 

(114) 

where K\ 2 — K\ + K 2 . By using Feynman parametrization we have 
h{ei) = / 4 [(ei,a)(e 2 ,a)(e 3 , a)(e 4 , a)] - ^(e^, e j )/f +2 [(e fc , a)(e i? a)] 

+ ((e 1 ,e 2 )(e3,e 4 ) + (ei, e 3 )(e 2 , e 4 ) + (ei, e 4 )(e 2 , e 3 )) / 4 D+4 [1]. (115) 
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We also assume that ei )3 satisfy the physical conditions in the two-mass- 
easy box integral and ei ;2 satisfy the physical conditions for the 2-mass-hard 
box integral. As we said before this can always be done by expanding all 
polarization vectors in terms of the 2 spinors from the 2 massless momenta. 

The direct calculation of the rational part of the two-mass-easy box inte- 
gral gives rather complicated formulas for generic polarization vectors (even 
after using the physical condition for €13). By choosing appropriate reference 
momenta for ei ;3 , the result simplifies greatly. In particular the reference mo- 
mentum of ei should be k 3 and the reference momentum of 63 should be k\. 
In some sense this is equivalent to the tensor reduction with the two factors 
(ei,p)(e 3 ,p — k\ — K 2 ). The explicit results are given as follows: 



j- / , r , r 1 [K% + s Kl + t\ , . w , s 

i4(A 3 Ai,e 2 , AiA3,e 4 ) = -- I R2 _ - + R2 _ - ) (e 2 , fci)(e 4 , h) 



1 (Kl + t K 2 A +s\ , , w , N 5 



4 \Kl-t Kl-sj ^^>^^> 9 



£2, fc 3 )(e 4 , h) - 77 (fa, k 3 )(e 2 , e 4 ) 



9 

~ 5 
^4(AiA 3 ,e 2 , AiA 3 ,e 4 ) = - (l|e 2 |3] (l|e 4 |3] 



+ I ((e 2 ,fa)(€ 4 ,k 3 ) + (e2,fc 3 )(e4, fci)), (116) 



9 

(1|X 2 |3] 2 



(e 2 , fa) (e 4 , fci) (e 2 , fc 3 ) (e 4 , fc; 



(X 2 2 - a ) (tf* - (K 2 2 -t)(K 4 2 - S ) 



(117) 



Other cases can be either obtained by conjugation or by relabelling fa 3 . 
In fact J 4 (AiA 3 , e 2 , A3A1, e 4 ) = / 4 (A 3 Ai, e 2 , AiA 3 , e 4 ) as it is invariant under 
conjugation. 

For the general case, we use the reduction formula by changing the refer- 
ence appropriately. For example we have: 

^4(^1,62,773X3,64) = I 4 {fa,e 2 ,e 3 ,e 4 ) + -Jj^ / 4 (ei, e 2 , £; 3 , e 4 ) 

7 4 (A3Ai, e 2 , A1A3, e 4 ) 



(13) 



2 



+ iVl ^^ 1] h(h,e 2 ,h,e 4 ). (118) 



34 



By using the explicit result of the 2-mass triangle integral we have 
K 2 + s 



I 4(^1, e 2 , ^3, e 4 ) 
+ 
+ 
+ 
+ 



K 2 A +t 



6(K% - s) 2 ^ U K2 ^ e2 > /Cl ^ 64 ' kl "> 
ei,K 4 )(e2,h)(e4,h) 



6(K 2 - 1) 2 

1 (K 2 + s K 2 + t\ 

12 [ k 2 - s Kl-t ) (( ei ' e2 K e4 ' fcl ) + ( ei ' e4 )v e2 ' fcl ) 



6(K 2 - s 
(ei,K 2 



+ 

— s 

(62, K3j(e4, Kl) 



(e2,Ki)(e 4 ,K 3 ) + 



(£1,^4) 



6 (X| - s) 6 (fTf - t) 
(ei,e 4 ) 

( ei ,X 4 ) (ei,X 2 



Q(K 2 - 1) 

((C2, fci) (C4, #4) - (C2, #2) (e 4 , fci)) 



+ g((ei,C2)(e4,fc3) + (ei,e 4 )(e2,^3) + (e 2 , e 4 )(ei, fc 3 )) 



+ - (e 2 , fci) (e 4 , ^4 



Iff - t Iff - 8 



(119) 



Except for the last term, this formula is invariant under the interchange 
2 <-> 4 (s <-> i). By setting ei = fci we get 

J 4 (fci, 62, A; 3 , e 4 ) = ^(2(fci, fc 3 )(e 2 , e 4 ) - ((e 2 , fci)(e 2 , fc 3 ) + (e 2 , A; 3 )(e 2 , £4))). 

(120) 

This agrees with the result by direct computation by first doing the tensor 
reduction and then using the result for bubble integrals. These formulas are 
quite useful for obtaining compact analytic formulas for QCD amplitudes. 

For easy reference we also give here the relevant formulas in terms of 
Feynman parameters: 



f 4 [(ei,a)(e 2 ,a)(e 3 ,a)(e 4 ,a)] 

(h,K 2 ) {h,K 2 ) ( (e 2 , h) (e 4 , ki) 



+ 



/f +2 [(e 2 ,a)(e 4 ,a)] 



(C2, fa) (€4, fa) 

SK 2 -s)(K 2 -t)^ {K 2 -t){K 2 -s) 
1 



6{K 2 + K 2 -s-t) 
(e 2 ,/ci) (e 4 , fci 



(e 2 , fcx) (e 4 , fc 3 ) + (e 2 , fc 3 ) (e 4 , fcj) 



^2^2 1 ' ' |M - ! ; 1^2, fc 3 ) (64, fc 3 J \1 , v 

( 2 4 " st) v ( ^- s) (Ki-t) + (K|-t)(^- s )JJ- (121) 
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All other I 4 +2 [(cj, a)(e,-, a)]'s are identically zero. In the above formulas e 13 
satisfy 2 conditions: (q, fci j3 ) = 0, i.e. the reference momentum of l\ is A; 3 
and the reference momentum of e 3 is k\. We also note when (ei,e 3 ) is not 
equal to 0, it cancels with the factor (K 2 + K\ — s — t) = (ki, k 3 ) = — (ei, e 3 ). 

7.4 The box integrals of degree 4 polynomials: the 
two-mass-hard case 

As before, the direct computation of the rational part of the two-mass-hard 
box gives very complicated algebraic expressions. In fact it gives the most 
complicated formula up to now. Changing the reference momenta simpli- 
fies a little bit, but the resulting formula is still not intelligible. The main 
complication comes from the presence of the 3-mass triangle integrals. To 
organize the final result, we proceed to do tensor reduction one more time. 
To begin with let us define the integral we want to compute: 

ll mh {ei,e 2 ,e 3 ,e 4 ;c 3 ,c A ) 

= h[(e 1 ,p)(e 2 ,p - k 1 )((e 3 ,p + K 4 ) + c 3 )((e 4 ,p + K 4 ) + c 4 )] 

f d D p (e l7 p)(e 2 ,p- k l )((e 3 ,p + K 4 ) + c 3 )((e 4 ,p + K 4 ) + c 4 ) 
~ J m D / 2 P 2 (p-k i y(p-k 12 ) 2 (p + K 4 ) 2 

(122) 

The external momenta are k\, k 2 , K 3 and K 4 as shown in Fig. 8 by setting 
K\ = k\ and K 2 = k 2 . k± t2 are the two massless external legs and K 3j4 
are the two massive external legs. We set t — (k 2 + K 3 ) 2 = (K 4 + ki) 2 . 
We also require that e^ 2 satisfy the physical condition, i.e. (e±,ki) = and 
(e 2 , k 2 ) = 0. The formula for generic 6i >2 is beyond the scope of this paper. 
(e 3j4 are arbitrary 4-dimensional polarization vectors.) 

There are 4 possible cases for e\ y2 . In the same helicity cases, the box 
integrals can be easily reduced to triangle integrals by using the reduction 
formulas eq. (14) and (16) in Sect. 3. So we will only consider the difficult 
cases where ei >2 have different helicities. To be definite we set t\ = Xifji 
and e 2 = i] 2 X 2 . The opposite case can be obtained from this one simply by 
conjugation. We will give the explicit formula for this case at the end of this 
subsection. 

To do tensor reduction for the factors associated with the two massless 
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external legs, we have: 

T + = (ei,p)(e 2 ,p-fci) = <2|AT 3 |1] ' ^ ^ 

The middle factor in the above can be decomposed by moving the first factor 
of (p — k\) towards the second (p — fci). Explicitly we have: 

(p - h) k 2 K 3 h (p - h) —T± + T 2 , (124) 

7\ = I^( P -k 1 )k2K 3 + I^(K 3 k 1 -k2K 3 )p + I^K A k i p, (125) 
T 2 = pL^Ksh + il^ -t)k 2 k lP 

+ I^{-I^K 3 + I^\k 2 + K 3 )~(I^ -t)k 2 ). (126) 

Omitting the factor (2|fT 3 |l], we now compute the various terms. We have: 

{V2\T 2 \fn) = p^ri^hK-shlf)^ 

+ (J« - t)((rj 2 2) [fk l](AiA 2 ,p) - / (2) faM^i)) 

+ /( 1 )(/( 3 )(r ?2 |(A; 2 + K3)|^)-/ (2) fa|ir3|^i)). (127) 

This gives the following rational terms: 

A 2 = -^(v2\k 2 K 3 k 1 \f ll )(e 3 ,e 4 )-t(r ]2 2)[fi 1 l]li m \\ 1 \ 2 ,e 3 ,e 4 ) 

+ t ( V2 \k 2 \^] J 3 3m (e 3 , e 4 ) + ( V2 2) [77! 1] Q«l|e 3 |2] c 4 + (l|e 4 |2] c 3 ) 
+-L«l|e 3 |2] e 4 + (l|e 4 |2] e 3 , 7 fci + 2 fc 2 + 9 tf 4 ) 

lo 

+ <7/2 2) [771 2] ((e 3 , k 12 ) (e 4 , fc i2 ) - 2 s 12 (e 3 , e 4 )) 
lo 

+ (?7 2 |(fc 2 + K 3 )|?7i]((e 3 , fc 2 + X 3 )(e 4 , k 2 + K 3 )-2t (e 3 , e 4 )) 

- ^(r ?2 |K 3 |r7 1 ]((e 3 ,K 3 )(e 4 ,K 3 )-2^3 2 (e 3 ,e 4 )) (128) 

The explicit formulas for the other 3 terms in T\ are not quite illuminating 
and we refrain from writing the explicit results here. We can write the 
result in terms of the rational parts for the three-mass and two-mass triangle 
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integrals arising from the tensor reduction. We have: 

M = (2|^3|77i)(/ 3 2m (^A 2 ,e 3 ,e 4 ) 

+ / 3 2m (^A 2 , (c 3 - (e 3 , K 3 ))e, + (c 4 + (e 4 , ^4 + fci))e 3 )) 

+ /3 3m (u,e3,e 4 ) + / 3 3m (u,(c a -(e3,iir3))e4 + C4e3), (129) 

+ (l fe |/f 4 |l](i? m (Ai^i,C3,C4) 

+ J 3 2m (A 1 r7 1 ,(c 3 -( e3 ,A; 2 + ir 3 ))e 4 + (c 4 + (e 4 ,ir 4 ))e 3 ), (130) 

where 

u = <7 72 j J ftr 3 |l]A 1 77 1 + falKs^Xw - (k 2 , Ks)^. (131) 
Combining the above results together we have: 

ll mh (Xifji, r] 2 X 2 , e 3 , e 4 ; c 3 , c 4 ) = ^r±^y. (132) 

(2|A 3 |1J 

For the opposite helicity case the formula is: 
ll mh (r]i\i, X 2 fj 2 , e 3 , e 4 ; c 3 , c 4 ) 



(^ilfei-fi^l^) (e 3 ,e 4 ) 

o 



(l|^s|2] 

- i (77! 1) [773 2] / 2m/i (A 2 A 1; e 3 , e 4 ) + * <^i|^ 2 |^ 2 ] /f m (e 3 , e 4 
+ (7 7l l>[77 2 2]Q((2|e 3 |l] C4 + (2|e 4 |l]c 3 ) 

+-L«2|e 3 |l] e 4 + (2|e 4 |l] e 3 , 7fci + 2A; 2 + 9KT 4 ) 

lo 

+ ^ 2 ) 2 ] (( e 3> ^12) (e 4 , few) - 2 S12 (e 3 , e 4 )) 
lo 

+ — <77i|(^2 + K 3 )\fj 2 ]((e 3 , k 2 + K 3 ){e 4 , k 2 + K 3 )-2t (e 3 , e 4 )) 

- 1 tf 3 |i»] ((e 3 , X 3 )(e 4 , tf 3 ) - 2K 3 2 (e 3) e 4 )) 
+ (77 1 |X 3 |2](7 3 2 -(A 2 772,e 3 ,e 4 ) 

+/ 3 2m (A 2 r/ 2 , (c 3 - (e 3 , ^ 3 ))e 4 + (C4 + (e 4 , K 4 + fci))e 3 )) 
+ /|™({), e3 , e4 ) + /3 3 -({),( C3 -( e3 ,K 3 )) e4 + C4e3 ), (133) 

+ (l|^ 4 |77 2 )(/ 3 2m (r /l A 1 ,6 3 ,6 4 ) 

+ I^imXi, (c 3 - (e 3 , A; 2 + K 3 ))e 4 + (c 4 + (e 4 , X 4 ))e 3 )l , (134) 



where 



v = {l\K 3 \rh) 771A1 + (2|X 3 |77 2 ) r]iA 2 - (fc 2 , i^r^. (135) 
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8 Extra terms for box and triangle tensor re- 
duction 



8.1 Extra terms for box tensor reduction 

In order to get a comparatively compact analytic formula for the two-mass- 
hard box integral, it is necessary to do tensor reduction because a direct 
computation of the rational part by the recursive method gives a very com- 
plicated formula. It would be a better idea to organize it into lower degree 
and lower point integrals. A naive tensor reduction directly in D = 4 would 
give an incorrect result because the box integral is ultra-violet divergent for 
a degree 4 polynomial of momenta in the numerator. It turns out that the 
difference between the D = 4 tensor reduction and the correct tensor reduc- 
tion is just a rational function. In this subsection we will compute this extra 
rational function explicitly. 

For a degree 2 polynomial g(e,k,p) in the internal momentum p, the 
general form of the D = 4 tensor reduction we used is as follows: 



We assume that all polarization vectors e and momenta k are 4-dimensional 
and the above relation is derived by assuming p is also a 4-dimensional mo- 
mentum. We also assume that g and g depend on the momentum p through 
scalar products (e,p) and/or (k,p). Because we use FDH regularization [72], 
p is actually promoted to be in D = 4 — 2e dimensions in our later calcula- 
tions of the rational part. This affects only the last term in eq. (136). For 
arbitrary dimensional internal momentum p, the above formula is still valid 
if we make the substitution: 



That is, the extra dimensional part (D — 4 = — 2e) of the momentum p must 
be subtracted from p 2 which actually stands for the scalar product in D 
dimensions in the subsequent computation in the four-dimensional helicity 
regularization scheme. Because pentagon and hexagon diagrams are ultra- 
violet convergent we can safely discard this term in the tensor reduction for 
5- or 6-point diagrams. This term does give a non-vanishing contribution to 
the rational parts for box and triangle integrals. 



g(e,k,p) = g(e,k,p) -p 2 f(e,k). 



(136) 



2 2 2 

P ^P -P D -4- 



(137) 
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For box integral we have 

= / dV(l-£<0 / dTT 3 
Jo i Jo 

*j^P\M,P)^P)* p2T - Ta - S -\ (138) 

by transforming it into a Feynman parameter integral and omitting terms 
which are vanishing in the limit e — > 0. 

The integration over the momentum p can be done easily by splitting it 
into a D = 4 part and a (— 2e)-dimensional part. We have then 

A hox = e (-(e 3 , e 4 )) /f +4 [1] = (e 3 , e 4 ) + 0(e). (139) 

By using this result the correct formula for computing the rational part of 
the box tensor integral is: 

h[g(e,k,p)((e 3 ,p) + c 3 )((e 4 ,p) + c A ))\ = I 4 [(g(e,k,p) - p 2 f(e,k)) 

x((e 3 ,p) + c 3 )((e 3 ,p) + c 4 ))}-lf(e,k){e 3 ,e 4 ). (140) 

o 

8.2 Extra terms in triangle tensor reduction 

For triangle integrals, using the same reduction formula as given in eq. (136) 
also gives incorrect results. The extra terms are given by the following for- 
mula (p is the internal momentum between the external legs K\ and K 2 ): 



I 3 [g(e,K,p)((e 3 ,p) + c 3 )} = I 3 [(g(e, K,p) - p 2 f(e,k)) 



x((e 3 ,p) + c 3 )] - /(e, k) 



X - c 3 + i (e 3 , K x - K 3 ) 



(141) 



We note that the rational part I 3 [g(e, K,p) ((e 3 ,p) + c 3 )] is defined by the 
Feynman integral without the usual minus sign, just as we did before in 
eq. (52). 
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9 Conclusion 



The calculation of multi-leg one-loop amplitudes in QCD is highly motivated 
through the need of precise predictions for multi-particle scattering at TeV 
colliders like the Tevatron and the LHC. Whereas unitarity based methods 
are very successful to provide the cut-constructible parts of one-loop ampli- 
tudes, the evaluation of rational terms which are not defined by the cuts 
alone, is a much harder question. Recently the bootstrap recursive proposal 
of Bern, Dixon and Kosower [34] has been made for their evaluation. (It was 
used to do real QCD calcuations [34, 70].) In this paper we have developed 
a general formalism targeting only the computation of the rational terms by 
using Feynman diagrammatic methods. 

In comparing with the previous use of diagrammatic methods in the com- 
putation of multi-leg one-loop amplitudes, the computation of the rational 
terms only is greatly simplified by using the BDDK theorem [3]. By using 
this theorem and a thorough analysis of the recursive relations we derived 
recursive relations also for the rational terms. Quite explicit results are given 
for all bubble and triangle integrals, and box integrals up to two-mass cases. 
These results are all the ingredients for computing up to the 6-gluon one-loop 
amplitudes in QCD. We will compute the (rational parts of the) 5-gluon and 
6-gluon amplitudes in [68] and [69] respectively. 

By using the recursive relations it is straightforward to derive the rational 
terms of the 3-mass and 4-mass (tensor) box integrals. The formulas obtained 
are quite complicated. It may only be useful for implementing numerically. 

The method developed in this paper can be extended to the cases with 
massive internal particles and internal (massive) quarks. Our method com- 
plements the twistor-inspired approach quite nicely. We envisage that an 
automatic implementation of both approaches could push the present limit 
of 6-gluon to about 8-gluon or higher. 
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